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Abstract 

Multiple parton scatterings inside a large nucleus generally involve higher- 
twist nuclear parton matrix elements. The gluon bremsstrahlung induced by 
multiple scattering depends not only on direct parton matrix elements but 
also on momentum-crossed ones, due to the Landau-Pomeranchuk-Migdal in- 
terference effect. We show that both types of twist-four nuclear parton matrix 
elements can be factorized approximately into the product of twist-two nu- 
cleon matrix elements in the limit of extremely large nuclei, ^4 ^ oo, as 
assumed in previous studies. Due to the correlative nature of the twist-four 
matrix elements under consideration, it is actually the off-forward parton 
distributions that appear naturally in this decomposition, rather than the or- 
dinary diagonal distributions probed in deeply-inelastic scattering. However, 
we argue that the difference between these two distribution classes is small in 
certain kinematic regimes. In these regions, the twist-four nuclear parton ma- 
trix elements are evaluated numerically and compared to the factorized form 
for different nuclear sizes within a schematic model of the two-nucleon corre- 
lation function. The nuclear size dependence is found to be A^^^ in the limit 
of large A, as expected. We find that the factorization is reasonably good 
when the momentum fraction carried by the gluon field is moderate. The 
deviation can be more than a factor of 2, however, for small gluon momentum 
fractions, where the gluon distribution is very large. 
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Typeset using REVT^ 



I. INTRODUCTION 



The success of perturbative QCD (pQCD) in describing hard processes in high-energy 
colhsions is mainly attributed to the asymptotic freedom of QCD at short distances 
and to factorization theorems 0. Specifically, the cross sections of processes that involve 
large momentum transfers can be factorized into a convolution of perturbative hard scat- 
tering cross sections and nonperturbative parton distributions and fragmentation functions 
that contain long distance physics. Even though they are not calculable within pQCD, 
these parton distributions and fragmentation functions can be rigorously defined in QCD 
independently of any specific process and measured in many different experiments. Such 
factorization has been proven up to next-to-leading twist (twist-four) for hard processes 
involving both hadrons and nuclei. We will refer to this as the generalized factorization. 

The leading twist-four contributions in hard processes in nuclei normally involve multiple 
scattering with partons from different nucleons. They generally depend on twist-four nuclear 
parton matrix elements such as 

^ dy^dy^ e'^^^^y-+^^^^^^y^-y~^^e{-y^)e{y- - yf) 

x\{A\U^)^^F^+{y,)F+-{y-)Uy')\A) , (1) 

which describes the quark-gluon correlation in a nucleus. This matrix element also appears 
in both lepton-nucleus deeply-inelastic scattering (DIS) [|] and in Drell-Yan cross section of 
pA collisions [0,0. We will work in the infinite momentum frame, where the four-momentum 
of the virtual photon and the nucleus (atomic number A) have the form 

g= [-gV2g-,g-,oj, 

PA = A[p+,0,6j_], (2) 

respectively. The Bjorken variable is then xb = Q"^ /2p^q^ . Our convention for four-vectors 
is = [k~^, k~ , k±], where 

r.i^. (3, 

Assuming that the two gluon fields in the rescattering process associated with Eq. (|ip 
come from the same nucleon in the nucleus due to color confinement, it has been argued 
^ that the above twist-four nuclear matrix elements are enhanced by a factor of A^^^ as 
compared to the leading twist quark distributions in a nucleus, 

for A^ 1. For processes involving a large transverse momentum scale i"^ ^ Q"^, the ratio of 
the twist-four contribution and the leading twist one is therefore proportional to asA^^^ /i'^. 
For large values of A, where the above analysis is valid, this quantity can be related to an 
expansion parameter. In this sense, the above matrix element is the leading higher-twist 
contribution to hard processes involving multiple parton scattering in nuclei. 
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In a recent study ^ , Guo and Wang extended the generalized factorization approach to 
the problem of parton energy loss and modified quark fragmentation in DIS off a nuclear 
target due to gluon bremsstrahlung induced by secondary quark-gluon scatterings. Because 
of the Landau-Pomeranchuck-Midgal (LPM) interference effect, gluon bremsstrahlung 
with small transverse momentum, or large formation time (rj ~ Q'^/Mi'^ in the nucleus 
rest frame, where M is the nucleon mass), is suppressed. This limits the available phase 
space of the transverse momentum to i'^ ^ Q'^/MRa ~ Q'^ /A^^^, ensuring the validity of the 
leading logarithmic approximation in the study of jet fragmentation for i"^ <^ in a large 
nucleus with A^^^ ^ 1. The twist-four contribution to the modified fragmentation function 
in this case is proportional to agA^^^/Q"^, which depends quadratically on the nuclear size. 
Such a novel quadratic nuclear size dependence has recently been verified by the HERMES 
experiment |p!0[ . 



Similar to other twist-four processes in a nucleus, the nuclear modification to the frag- 
mentation function is also proportional to twist-four nuclear parton matrix elements. The 
LPM interference effect is explicitly embedded in the combined twist-four parton matrix 
elements. The quadratic nuclear size dependence of the modification to the fragmentation 
function is based on a generalized assumption that the twist-four parton matrix elements 
factorize into twist-two parton distributions in nucleons [^]. The same approximation has 
been assumed for the momentum-crossed twist-four parton elements of a nucleus. This is 
a crude assumption at best, and does not specify the condition of validity nor provide any 
insight into the the relationship between nuclear and nucleonic parton distributions, which 
should depend on the nucleon wavefunction inside a nucleus. 

It is difficult to determine the validity of this approximation within the framework of 
pQCD since nucleons do not appear explicitly in this theory. However, many hybrid models 
which employ nonrelativistic quantum mechanics to define a wavefunction for nucleons in 



a nucleus have been developed in the literature [fllHlt^. In these models, the wavefunction 
allows one to decompose nuclear parton distributions into nucleonic ones. Phenomenolog- 
ically, these models assert that scattering processes involving nuclei can be understood as 
weighted averages over scattering processes involving nucleons. Interactions among nucleons 



are reflected in the wave function. The simplest of these models arrives at the relation 

da 
a 



rA do, 

fa/A{x) = A — pia) fa/Nix/a) (5) 

J X OL 



between the leading-twist distribution of a-type partons in a nucleus of size A and the same 
distribution found in a nucleon. The correspondence is made via a light-cone nucleon density 
function which is the probability of finding a nucleon in a nucleus with longitudinal 

momentum fraction a, normalized to 1. Experimentally, Eq.(^) is approximately satisfied 
for X ^ 0.1. At smaller values of x, the phenomenon of shadowing prevents realization of 



this naive model 15.16 



For twist-four matrix elements, one expects similar results. As twist-four objects are 
associated with partonic correlations, one expects two types of contributions : effects associ- 
ated with partonic correlations within a single nucleon and those associated with nucleonic 
correlations within the nucleus. The former effects, which involve nucleonic twist-four dis- 
tributions, are a simple extension of the convolution model; one simply substitutes new dis- 
tributions for the twist-two ones. Contributions from nucleonic correlatons imply multiple 
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scattering within the nucleus, and are inherently new effects. In particular, the distribution 
p{a) in Eq.(|D will be replaced by a more complicated distribution describing two-nucleon 
correlations, or the momentum-sharing between nucleons in the nucleus. Since these latter 
effects involve two nucleons, they will be enhanced by factors of the nuclear radius relative 
to the nucleonic higher-twist effects. This makes them the dominant contributions in the 
limit of large A. 

The purpose of this paper is to show that an analysis of these effects in the spirit of 
the convolution model reveals contributions from a more general class of twist-two matrix 
elements. Specifically, the decomposition necessarily contains off-forward parton distribu- 
tions (OFPD's) [pT| , p^ rather than just the simple diagonal matrix elements. At present, 
deeply-virtual Compton scattering (DVCS) is the only known process which involves the 
OFPD's explicitly [jl8[. We will demonstrate that these elusive objects could in principle 



also be probed in multiple parton scattering processes in a nucleus. 

This paper is organized as follows. In Section II we will give a brief review of multiple 
parton scattering in nuclei and modified fragmentation functions, focusing on their relation 
to the twist-four parton matrix elements in nuclei. Section III presents the main result of 
this paper, in which we derive a relation between a twist-four nuclear matrix element and a 
convolution of two nucleonic OFPD's. In Section IV, we discuss the properties of our result 
in certain limits. We will find that an analytic relationship between the simple factorized 
expression mentioned above and our result is not obvious, but numerical models show that 
they produce the same results in the limit of sharply-peaked nuclear wave functions. Section 
V contains some discussion about and conclusions of our results. 



II. MULTIPLE PARTON SCATTERING 

DIS on a nuclear target is the simplest environment in which to study the problem of 
multiple parton scattering in a nucleus. In this quark is struck by an energetic 

virtual photon and then scatters again with partons from other nucleons inside the nucleus. 
The rescattering will induce gluon bremsstrahlung by the propagating quark and cause 
the leading quark to lose energy. Such radiative energy loss will be manifested in the 
modification of the quark fragmentation function as compared to the one measured in DIS 
off a nucleon target, where there is no such rescattering. The gluon bremsstrahlung will 
interfere destructively with the final-state radiation of the quark-photon scattering. This 
LPM interference effect will give rise to some novel nuclear effects in the modified quark 
fragmentation function. 

Applying the generalized factorization of twist-four processes to the exclusive process of 
hadron production in DIS on a nuclear target and performing a coUinear expansion with 
respect to initial parton transverse momentum, one can obtain an effective modified quark 
fragmentation function with leading higher-twist contributions 

Jo ZJC Jzh Z L 

+ Ajg^ggiZ, X, XL, iT)Dg^hiZh/ Z^ . (6) 

Here, Dg^h{zh, fi^) is the usual renormalized twist-two quark fragmentation function in 
vacuum that satifies the normal DGLAP [|l^] QCD evolution equation. The additional terms 
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FIG. 1. A central-cut diagram for quark-gluon rescattering processes. 



are the leading higher-twist contributions from multiple parton scattering and induced gluon 
bremsstrahlung. These contributions are very similar in form to the normal gluon radiation 
in vacuum except that the modified splitting functions, 



[l-z) 



1 + 



T^^^{x,xl) + 6{1 



^iNMx) 



(7) 
(8) 



depend on the twist-four two-parton correlation function 

dy- 



Tqg{x,XT,XL] 



■dy^dy,d{-y,)d{y- - y^){l - e"-P^(^r-J'"))(l - e^^^^^^^O 



where Xl = i'^/2p'^q~z{l — z) and Xt = (k^) /2p'^q^z. The virtual corrections supply the 
5-function contribution to the '+ '-function, along with the explicit end point contribution 



dz 



1 + z^ 



TggiX.Xl) 



dz- 



1-z 



2Tt{x,XL)l 



'19^ 



-il + z 



(10) 



required for conservation of quark flavor. 

The twist-four parton matrix elements are in principle not calculable and can only be 
measured in experiments, just like twist-two parton distributions. However, under certain 
assumptions, one can use some nuclear model to relate them to twist-two parton distributions 
in nucleons. Along the way, one obtains the A-dependence of these nuclear matrix elements. 

If we assume that the nuclear wave function can be expressed as a multiple-nucleon state, 
with each nucleon a color singlet, the two gluon fields must operate on the same nucleon 
state inside the nucleus. Consider the dominant case where the quark and gluon fields 
operate on different nucleons inside the nucleus. The integration over yjf and in Eq. (^ 
should give the length scale t^Ra, where r^r is the nucleon radius and Ra — 1.12^4^/^ fm 
is the radius of the nucleus. The twist-four nuclear parton matrix elements should then be 
approximately proportional to A^^^. If the quark and gluon come from the same nucleon, 
the matrix elements will only be proportional to A, which is subleading. 
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The twist-four two-parton correlation function that enters the modified quark fragmen- 
tation function has not only normal parton matrix elements representing direct terms in 
the square of the amplitude, but also those representing interference. The former has mo- 
menta flowing directly along two parton fields separately while the latter has momenta 
flowing across two different parton fields. These two different contributions were called 'di- 
agonal' and 'off-diagonal' matrix elements, respectively, in Ref. |^. We call them 'direct' 
and 'crossed' here to avoid confusion with truly off-diagonal matrix elements, in which the 
momenta of the external states are different. The relative signs between these two kinds of 
matrix elements reflect the physics of the LPM interference effect in the processes of induced 
gluon radiation. As illustrated by the central cut-diagram in Fig. |I], the gluon radiation can 
either be produced as final state radiation of the photon-quark hard scattering or initial 
state radiation of the quark-gluon rescattering. In the former case, the energy of the radi- 
ated gluon is provided by the initial quark with x = xb + xl- The quark-gluon rescattering 
can be very soft since the momentum fraction carried by the gluon, x^ = xt = (k^) /2p'^q~ z, 
is very small when (k^) ^ 0. This is normally referred to as a hard-soft process. In the 
latter case, however, the initial gluon must carry a finite momentum fraction X3 = xt + xl to 
induce the gluon radiation. Such a process is called double hard scattering. Contributions 
from central-cut diagrams such as Fig. Q contain both of these processes as well as their 
interference. Since initial and final state radiation amplitudes have a phase difference of 
TT , the final result of the sum of these contributions is the dipole-like form-factor in the 
radiation spectrum which is now absorbed into the definition of the two-parton correlation 
T^{x,xt,xl)- As xlP'^ — > 0, the effective gluon radiation spectrum vanishes because the 
interference becomes complete. This is exactly the LPM effect 0, which is now embedded 
in the effective two-parton correlation function. 

Because of the LPM interference effect in induced bremsstrahlung, the effective two- 
parton correlation function that enters the modified fragmentation function essentially con- 
tains four independent twist-four nuclear parton matrix elements. The two direct ones 
correspond to gluon radiation associated with photon-quark and quark-gluon scattering. In 
the corresponding forward scattering processes, shown in Fig. |l|, momentum flows separately 
along the quark and gluon lines. The twist-four parton matrix elements can then have the 
interpretation of a two-parton joint distribution inside the nucleus. The two crossed ma- 
trix elements are related to the inteference between the two different radiation processes. 
In this case, there is actually a momentum flow in the amount xlp'^ between quarks and 
gluons in the forward scattering amplitude. Such crossed matrix elements do not have the 
interpretation of a normal parton distribution. Since we only consider the case where the 
quark and gluon are from different nucleons inside the nucleus, these contributions should 
be related to the off-forward (or skewed) parton distrbution functions of a nucleon. In the 
limit of vanishing skewedness {xl — *■ 0), the crossed matrix elements approach the direct 
ones. One can see by inspection that the direct matrix elements are real, while the crossed 
ones are complex. 

In the rest of this paper, we study these twist-four nuclear parton matrix elements within 
the convolution model and relate them to generalized nucleonic parton distributions. The 
matrix element at issue, 

K{xi,X2,Xl) = / -;— dy^ dy^ e{-y:^)e{y' - y^) 
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FIG. 2. An example of a crossed parton matrix element. Note the momentum transfer from 
^luon to quark. 



describes the removal of a quark with momentum Xip and a gluon of momentum (x2 + Xl)p 
from our nuclear state \A), and the subsequent replacement of a gluon with momentum X2P 
and a quark of momentum (xi + as illustrated in Fig. 0. This nuclear parton matrix 

element is useful in constructing the physical combinations that appear in many nuclear 
scattering processes [Q, as well as the in- medium evolution of the parton fragmentation 
functions In particular, the correlation function (P) appearing in the medium-modified 
quark fragmentation function can be written as 

Tfg{x, Xt, Xl) = K{x + XL, XT, 0) - K{x, Xt, Xl) 

-K{x + xl,xt + xl,-xl) + K{x,xt + xl,0) . (12) 



III. THE CONVOLUTION MODEL 

In order to relate the matrix element ([TT|) to the nucleonic degrees of freedom of the 
nucleus, we must define these degrees of freedom quantitatively. Our formalism is based on 
light-cone perturbation theory [^], in which one expands the physical state under consider- 
ation in terms of the free-particle states of its constituents. In our case, we would like to 
express our nuclear state in terms of free nucleonic states. Phenomenologically, we know 
that this basis is not complete and therefore cannot be guaranteed to span our nuclear state 
space. 

Contributions from higher Fock states are required in quantum field theory to generate ef- 
fective nonlocal interactions from the underlying contact terms of the Lagrangian. However, 
most of the effect of higher Fock states is already included in the definition of the twist-two 
nucleonic parton matrix elements. The additional effects of higher Fock states in a nucleus 
are induced by nucleonic interactions. However, these effects are limited by the effective 
interaction energy of the system under consideration. Here, we are mainly concerned with 
nuclear systems whose interaction energy is small compared with the energy per nucleon. 
This allows us to consider only the lowest Fock state, that of A nucleons. We mention here 
that although we neglect interactions in the form of higher Fock states, nucleon correlations 
will still occur through the wave function in our Hilbert space. These correlations lead to 
nontrivial relations between the nuclear matrix element and the nucleonic distributions. 

In light of the above approximation, we write 
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\A) ^ J dUA mP^})\{P^}) 2p+(27r)3(5(3) (^PA - Ep^) ' (13) 
where represents the state of A free nucleons of momenta {pi} normahzed as 

mm) = Y[Mi^^)'^^'\p^-pi) , (14) 

i 

and 

represents the differential phase-space of A nucleons. 

The nucleon states are specified by their '+' and '_L' momentum components; d^pi = 
dpfd'^Pi±. The '-' component of each momentum is determined by the on-shell condition, 

_ pL + 

P. = , (16) 

and is not conserved (i.e. J2iP7 7^ Pa)- Again, M is the nucleon mass. The normalization 
of our nuclear state, = 2p^(27rY5^^\pA ~ p'a) iniphes 

/ dllA |0(fe})p2p+(27r)35(^) [pa - = 1 . (17) 

If we were to consider higher Fock states as well, the 1 on the right-hand side of this equation 
would be replaced by the probability of finding our state in this lowest Fock state. 

Our wave function contains many distribution functions. In particular, we can define 
the one-nucleon density 



pi{k) = J dUA-i |0(A;, fe})p 

x2p+{27rf6(Ap+ -k+ -J2pi]6'^'Uk^+Y.p,^] , (18) 



|2 



i=l / \ i=l 



which represents the probability of finding a nucleon of momentum k in our nucleus irre- 
spective of the momenta of the other nucleons. In light of Eq. (plTl), pi satisfies 

d^k 9(k+) 

^ ' Pi{k) = 1 . (19) 



(27r)3 2k 

Hence the light-cone nucleon distribution function in Eq.(|^) can be written as 

where a = k~^ /p'^. Eq.(^ comes directly from the substitution of ([T3|) into the definition of 
the nuclear quark distribution function fg/A^x) in Eq. (^). 
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The two-nucleon correlator, 



P2(A;i; A)= J dIlA-2 (l>*ih - A/2, h + A/2, {p,})(t>ih + A/2, h - A/2, fe}) 

X 2p^{2nf5 (^Ap+ - kf - k+ - Y.jt^ ^^'^ (^i^ + + E P^^) > (21) 

contains information about the sharing of momentum by nucleons inside the nucleus and 
appears in the double scattering process we consider. The two-nucleon density, ^2(^1; ^2; 0), 
represents the probability of finding two nucleons with the specified momenta within the 
nucleus, and is normalized as 

" dU2P2{k,;k2;0) = l . (22) 

The use of this function rather than the two-parton correlator, p2, leads to the expectation 
of diagonal parton distributions in twist-four nuclear matrix elements. 

Due to the ^-functions, our matrix element cannot readily be interpreted as a product of 
twist-two distributions as it stands. These 6'-functions order the fields along the light-cone 
axis to make the multiple scattering process physical. Employing the representation 

for the 6'-function, our nonperturbative distribution takes the form 

1 f dzi dz2 

zi — uji — ie Z2 — UJ2 + is 

rlii rill, rlii„ ^ip+y- (^^^j^^i-uji) ^ip+ {x2+xl-zi+uji) ^ip+y~ {z2-uj2~X2) ^24) 



K{xi,X2,Xl) = - J 



2tx 27r 27r 



A 



Here, uJi and ljJ2 are arbitrary real variables that will be chosen later to simplify the final 
results. 

Substitution of our approximate nuclear state (|l^) into ( p^ ) leads to matrix elements of 
the form 

m\ml^i^{y')K^{y2)K''{yl)\{p^}) ■ (25) 

Assuming that the color correlation length along the light-cone within our nucleus is not 
larger than the nucleon size and neglecting the effects of direct multi-nucleonic correlations 
as higher twist, we can factorize this expression into a product of single-particle Hilbert 
space amplitudes: 

A{p[\mi^^{y-) F;+(y2-)Fr(l/r)bi) \{2pt{2nf5\p[-p,) 

i=2 

+ (^^{p\\mi^i^{y-)\pi) {P2\F:Hy2)F:-iyi)\P2) l[2pt{2nr6%p[-p,) . (26) 



9 



Our dismissal of direct multi-nucleonic correlations is one of the main approximations in 
this paper, and will be discussed in more detail in the conclusions. For now, we note that 
these higher-twist corrections are suppressed by powers of and as such can be neglected 
at large scales. 

The contributions to (|26D can be interpreted in terms of the multiple parton scatter- 
ing picture in DIS. Matrix elements related to '?/'(0)7"^^/'(?/~) represent the probability that a 
quark in a certain nucleon is struck by our probe. The struck quark then propagates through 
the nucleus, encountering another parton at some point during its journey. If this rescat- 
tering occurs while the struck parton is still in its parent nucleon, the effect is represented 
by a twist-four nucleonic matrix element convoluted with the single-nucleon density. This 
is essentially the first term in Eq. fpBI) . Considering the coherence length of the scattering 
to be of the order of the nucleon size, this term should be proportional to A. 

If the rescattering is with a parton in another nucleon, the probability is related to the 
twist-two parton distributions in each nucleon convoluted with the two-nucleon correlator. 
We can call this double-factorized rescattering. Since the nuclear radius grows with A^/^, 
one expects this double scattering process to be proportional to A^l"^ , enhanced by A^l^ over 
the double scattering within the same nucleon. The first term in the above decomposition 
is therefore suppressed relative to the second as A increases. 

Using the above decomposition, and keeping only the contribution of the double- 
factorized rescattering, K{x\,X2-,xj)) becomes 



i^(Xi,X2,XL 




(fki^ (fk2± <f^± 
zi-uji- le Z2-uJ2 + ie J {2tiY {2Txf (27r)2 
ai + C) e{ai - C) e{a2 + C) d{a2 - C) 



277 27r 27r 2(ai + C) 2(ai - C) 2(^2 + C) 2(a2 - C) 
X 5{xl - 2C + Z2 - zi - U2 + uji)p2{ki; k2; A) 



(27) 



X 



dy 

2n 



Jp+y (xi+zi-oJi+C) 




27C p^ 



Jp+y, {x2+XL-zi+uJi-(;) 



k2 + 



A 



Vd 




where = yi - y2, ai = kf/p+ and C = -A+/2p+. 

If we could take A = 0, this expression would be reduced to the ordinary parton distri- 
butions probed in deeply-inelastic scattering. As it is, we can employ the off-forward parton 
distributions (OFPD's) [jl^ [Note that the distributions used here are slightly different than 
those defined in [|l7|. In particular, G{x,^,t) = xFg{x,^,t)]. Thus, 



dy 
2-n 



oW y ^ 



ki 



dy 

271 

where 



k2 + — 



A 

pa 



y~ 



7+^ 




2Q (^., 



My 



(28) 
(29) 



(30) 
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represents the squared four-momentum transfers for the two matrix elements. Note ti t2 
since the '-'-components of our four- vectors are not conserved in this version of perturbation 
theory. 

In terms of these OFPD's, we have 



K{xi,X2,Xl) 




dai da2 d( 6{ai + Q 6{ai — Q 6{a2 + C) ^(ct2 — C) 



^P2ih;k2;A) (31) 



(27r)2 (27r)2 (27r)2 

7 ~ ^ 7 , . , S{Xl -2C + Z2- Zi~UJ2+ UJi) 

[Zi -UJx- le) [Z2 -UJ2 + i£> 

fxi + Zi-UJi + ( ti\ fx2+XL-Zi+UJi-C t2 \ 

^"^^i a, '-^'M^j^i a2 '^'M^j " 

This comphcated expression can be reduced to a more enhghtening form through a few 
plausible assumptions on the form of our nucleonic correlation p2 and the OFPD's. To begin 
with, we assume that p2 is peaked around = 1 and ( = A± = ki± = 0, with widths that 
are governed by the nuclear radius, R^. This ansatz is dictated by the expectation that the 
nucleons are confined within the nuclear radius in position space. Specifically, we write 




P2(ai,/ci±;a2,/c2±; -2C, A_l) = ra 

\ 

(32) 

with xa = ^/{MRa) and r2 approximately independent of A. The behavior of the normal- 
ization R^/x'a as a function of A can be determined from the normalization condition (p2D 
on our two-particle density. Using the (5-function to perform the ( integration and making 
the changes 

(33) 
(34) 
(35) 
(36) 
(37) 

in variables, we arrive at the expression 

I fA\ r . . 1 1 





= (oi - 1)/xa 


Vil. 


= RAki± 


s± 


= RaA± 


z 


= {zi + Z2)/2xa 


u 


= (zi - Z2)/2xa 



K{xi,X2,Xl) 



AttR\ \2 J J z + u + ^ — is z — u — ^ + is 

dui du2 6{1 + xa{j-'i — u)) 0{1 + xa{j-'i + u)) 



27r 27r l + x^(z/i— -u) 1 + Xa{i^i+u) 

e{l + XA{v2-u)) e{l + XA{v2 + u)) 

1 + xa{i'2-u) 1 + xa{i'2 + u) 
J 7^7^7^^2(^^i,t^i±;^^2,t^2±;2M,5x)(l + XAZ^2) (38) 



(27r)2 (27r)2 (27r) 
+ XAi^ + z) 

1 + xai^i ' -^--'--j- \ 1 + XaZ/2 



/Xi+X^(^ + Z) A ^ fx2+ XAi^ - Z) ~ 

X Q — ,xau, ti\G \ , -xau, t2 
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for (|3T|). For simplicity, we have chosen ui = —uj2 = — xl/2 and defined the parameters 
^ = xl/C^xa) and 



ti 



(1 + X^Z^j)^ — x^u"^ 



Au^ + ((1 + XAl^i)^!. — '2XAUVi±y 



(39) 



Working from this form of K{xi, X2, xl), it is easy to see the large-A enhancement of the 
multiple-scattering contribution. Since r2 depends only weakly on A, this contribution to 
K scales like A^^^ as expected. We can simplify our expression further by assuming that r2 
is sharply peaked in the sense that all of its moments are finite. In this case, we can expand 
our integrand about the peak in r2 in the formal limit A —>■ oo and drop all non-leading 
terms : 



K{xi,X2,Xl) 



A' 



S7iR\ 



du dz- 



1 



1 



z + u + — is z — u — ^ + is 
X Q{xi + xaH + z),xau, 0)G(x2 + xaH - z), -xau, 0) 



where 



r2[u) 



dui di'2 d'^vi± d^V2L d'^5± , ^ r \ 
2 77r-T^^2(z^i, z^2, f2±; 2m, 6±) 



2tx 2ti (27r)2 (27r)2 (27r)^ 



(40) 



(41) 



Expression ( ^OD is the main result of this paper. Its derivation requires only the assump- 
tions that the lowest Fock state dominates the nuclear wave function and that the nucleonic 
correlator is sharply peaked with a width dictated by the nuclear radius. Strictly speaking, 
this expression is valid only in the formal limit A — > cxo. For any finite value of A, one must 
investigate the size of the derivatives of the OFPD's in relation to A^f^. While this investiga- 
tion must be done within a specific model, the contributions are expected to be small as long 
as the singular regions are avoided. Our implicit assumption that the OFPD's are analytic 
functions of the virtuality of the momentum transfer, tj, is supported by studies of these 
functions []T7ip0i . On the other hand, the kernels dictating the evolution of these functions 
(cf [0) imply that the OFPD's are not analytic functions of their second argument. This 
is why we have not expanded our integrand about xau = 0. 

To put this expression into a form suitable for numerical evaluation, we write 



1 



1 



z + u + ^ — ie z — u — ^ + ie 



1 

Yz 



1 



+ 



1 



z + u + ^ — ie z — u — ^ + is 



(42) 



and explicitly separate the w-integration into its pole and principal value parts. After a 
change of variables, this leads to 



K{xi, X2,Xl) 



A' 



f"^ dz r°° dw 



87iR\ 



[F{z,w-z-0-F{z,w + z + 
lo 2z Jo w 

+F{z, w-z + 0- F{z, w + z-^)] 

/■°° dz 1 
-[F{z,z + 0-F{z,z-0]j , 



(43) 



where 
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F{z, u) = [Q{xi + + xaz, xau, 0)G{x2 + — xaz, xau, 0) 

+Q{Xi+ Xl/2- XaZ,XaU,0)G{x2 + Xl/2 + XAZ,XAU,0)]f2{u) . (44) 

In deriving this result, we have used the facts that r2{u) = r2{—u) and that the OFPD's are 
even functions of their second argument, as can be seen by inspection. 

Several features of this expression are worth pointing out. First, we note that the imag- 
inary part of our matrix element is odd in ^. This causes the combination in Eq.(|l^) to be 
real, as required. One can check that the two-parton correlation T^( expressed 
in Eq. ([T^), depends only on the real part of K{xi,X2,xl). In addition, this combination is 
necessarily positive. The fact that these consistency requirements are satisfied is gratifying, 
but not unexpected. On the other hand, the relationship between our expression and the 
naive expectation 

K{xi,X2,0) r^Q{xi)G{x2) , (45) 

where Q{x) = Q(xi,0,0) and G{x) = G{x,0,0) are the ordinary diagonal distributions, is 
entirely unclear. Ignoring, for the moment, the non-analytic nature of the OFPD's, one 
can imagine expanding each term in the integrand about the peak of the nuclear sharing 
function r2 and dropping all higher-order terms. This reduces the OFPD's to ordinary parton 
distribution functions, bringing us closer to (|45|) . However, we cannot formally reproduce 
this simple dependence because of a remaining convolution between the distributions. This 
convolution causes our matrix element to sample the entire parton distribution functions, 
regardless of the values of Xi, X2, and xl- Its presence is a direct consequence of the 
correlative nature of the matrix element. Mathematically, it comes from the ^-functions. 

We could attempt to remove the correlation by expanding our result about xa = 0. 
However, corrections to the leading term in our expansion diverge, indicating nonanalyticity 
at Xa = 0. This does not mean that the leading term is not a good approximation to 
the full solution, only that the corrections cannot be expressed in terms of powers of xa- 
Nevertheless, assuming that the corrections are small for some range of xa near zero, one 
can compute Eq. (^) via contour integration directly in the limit of a;^ — * 0, but with fixed 
value of ^ = xl/^xa- The result. 



K{xi,X2,Xl) 



XA^O 



Ko{xi,X2,Xl) 

Q(si + xl/2)G{x2 + Xl/2) 
i /■°° du 



X m) + - T- [hiu +0- h{u - 0] , (46) 

( vr JO M J 

is of the form of the naive expectation. In order to determine the validity of this approxima- 
tion, we must calculate K numerically in some model and compare the results to the above 
approximation. 



IV. NUMERICAL RESULTS 



In this section, we explore some of the properties of expression ( ^OD in a specific model. 
While this model is certainly not expected to conform to the quantitative details of the 
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realistic nuclear wave functions, we expect its general features to be echoed in more realistic 
treatments. 

For simplicity, we assume a Gaussian form for f2 : 



= ^^e-"' . (47) 



The normalization is determined via (^) in conjunction with (|32D by assuming a Gaussian 
dependence of r2 on 6±. The constant (R^A^) is a measure of the transverse momentum- 
sharing among nucleons in our nucleus, and is expected to be of order one. Expanding 
our integrands about the peak in f2 and performing the integral over u, we arrive at the 
expression 

X I [D{z + + D{z - e)] + !^ [e-(^+«' - e-(^-«)'l | . (48) 



Here, 



f{z, 0) = Q{xi + xl/2 + xaz)G{x2 + xl/2 - xaz) 
+Q{xi + Xl/2 - xaz)G{x2 + xl/2 + xaz) 



(49) 



and 

D{x) = r rfte*'-"' (50) 



is Dawson's integral. Since D{x) — l/2x as x ^ oo, K{xi,X2,xl) is obviously non-analytic 
in xa at xa = 0, as mentioned above. However, it can be checked explicitly that this 
expression reduces to ( ^BD when Xa ^ 0. 

In deriving Eq. (^), we have assumed xau ^ 1 in the expansion around the peak of 
r2{u). Such an approximation does not necessarily represent the leading behavior of the 
integral. The pole contributions occur in the region m ~ ^, where xau ^ xl/2 can in fact 
be of the same order or larger than the first argument of our OFPD's, 

Q{xi + Xl/2 ± xaz, Xl/2, 0)G{x2 + xl/2 =f xaz, xl/2, 0). (51) 

In this case, the full generalized distributions are sampled. However, in the relevant region 
of the 2;-integration in Eq. (^) with a Gaussian form of r2{u), \z\ < ^. This causes the 
first variables in the OFPD's to be bounded by xi and X2, respectively. According to 
model studies of the OFPD's |]20|j2T| , the OFPD's can be approximated by the ordinary 



parton distributions when the first argument is larger in magnitude than the second. In 
addition, since the OFPD's are continuous and G{x,^,t) is expected to be positive definite, 
the variation of G{x,^,t) is small in the region x < |^| <C 1. Hence we expect very small 
deviations from the ordinary gluon momentum distribution even in the region X2 <iC Xl <^ ^■ 
The singular nature of the quark distribution function Q{x) at x = leads to large variations 
of the associated OFPD in the region xi < xl- This will essentially be the limitation of 
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FIG. 3. The dependence of the ratio of K to Kq on shown for the nuclei '^^S, ^^Ni, and ^''^Pb. 
The dashed hnes show the saturation ratio for each nucleus (xi, — > 0). We note that although the 
saturation ratios are not very close to 1, the curves are quite flat when xa > xl, and the ratio 
increases with nuclear size. 
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our approximation here. Outside of our applicability region, the full generalized parton 
distributions are needed to predict the behavior of the matrix element. 



Using the CTEQ parameterization of parton distributions from data p2[ , we can calculate 



our expression numerically and see how well the results follow the factorized form Kq. Since 
only the real parts of K enters the two-parton correlation T^, we will concentrate on the 
real part of our matrix element. The nature of (|^) is such that our matrix element samples 
the parton momentum distribution at all values of x rather than just those close to xi^2- 
Since the parton distributions are not known for x — > 0, we assume a simple extrapolation 
with a constant value for gluons and set the quark distributions to zero beyond the region 
of parametrization. The errors introduced by such extrapolation are negligible since the 
contributions to the integral from this region is very small. 

Figure ^ shows the ratio K{0.3, 0.01, xl)/Kq{0.3, 0.01, xl) as a function of xa/xl for 
three different values of xa- Since xl measures the momentum sharing between nucleons 
in our matrix element, it should be smaller than the characteristic momentum fraction, xa, 
in our nucleus. If xl becomes too large, the momentum transfer is suppressed beyond the 
simple exponential suppression in Kq. This explains the behavior of the curve for small 
xa/xl- P^s xa/xl increases, the ratio approaches the direct contribution represented by 
the dashed line in the figure. As can easily be seen from the graph, there is some residual 
dependence on the nuclear size. However, the saturation ratio changes by less than 15% as 
one changes A from 208 to 32. 

To study the residual dependence of our matrix element on the nuclear size, we plot 
the ratio K/Kq at ojl = as a function of xa for three different values of xi in Figure ^. 
The ratio is seen to drop monotonically as xa is increased, but the close proximity of the 
curves Xi = 0.2 and Xi = 0.3 indicates only slight dependence on xi when xi is moderate. 
As Xl becomes smaller, the dependence on both xi and xa becomes far more dramatic. 
In either case, the xa dependence is approximately linear. Since we have already dropped 
many terms of this order in xa, this behavior is not at all surprising. However, it can lead 
to large corrections to the factorized form for real nuclei (where xa is of order 0.04). 

The Xl dependence is clearly illustrated in Figure ^ Here, we see that the factorized form 
Kq is indeed a good approximation for moderate Xi. In the region 0.15 ^ a;i ^ 0.75, the ratio 
changes by approximately 15% for ^°^Pb, while in the more restricted region 0.2 ^ si < 0.7, 
the change is less than 8%. For smaller nuclei, the region of 'moderate xi is more restricted. 
The behavior of our curves for small xi is due to the fact that the denominator of our ratio 
diverges as Xi — >■ 0, while the numerator samples the distribution function and smears the 
divergence. For large Xi, the behavior is quite similar to that of the ordinary convolution 
model, Eq.(^), which is due to the fact that the denominator vanishes as Xi — 1 while 
the numerator remains finite because the partons share momentum. This effect is usually 
referred to as Fermi motion. The dependence of the ratio on nuclear size is not too large 
(on the order of 15% for moderate Xi) in the range 0.032 < < 0.059 considered. The 
tremendous dependence of the placement of these curves on the value of X2 can be attributed 
to the turnover of the gluon momentum distribution at small x. Since this turnover occurs 
at small x, we expect large variation only when X2 is extremely small. 

To see this explicitly, we plot the ratio K/ Kq as a function of X2 in Figure |[ Here, we 
can clearly see that the large variations are confined to X2 ^ 0.01. Changing the scale of 
the input distributions allows us to explore the behavior of our ratio as a function of Q^. 
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Ratio vs X 
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FIG. 4. K/Kq versus xa for three different values of xi. The xa dependence is approximately 
linear for small xa, with a slope whose magnitude decreases as xi increases. Although this de- 
pendence is nominally of order xa, it can lead to quite large corrections for real nuclei if xi is too 
small. 
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Ratio vs X 
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FIG. 5. K/Kq versus xi for ^^S, ^^Ni, and ^^^Pb. Although there is strong dependence, the 

ratio is quite flat for moderate xi. The strong dependence of the placement of these curves on the 
value of X2 is attributed to the decrease of G{x) for extremely small values of x. 
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10""* lO""* 10^ 10"' 

^2 

FIG. 6. The dependence of our ratio on X2 for two different values of . The increase in 
the ratio as X2 decreases at 2.5 GeV^ is attributed to the turnover of the input gluon momentum 
distribution. Since this turnover is pushed back beyond our cut-off at X2 = 10"^ for = 5.0 
GeV^, the increase is not present at this scale. In both cases, the dependence of the ratio on X2 is 
quite moderate for X2 > 10^^. 



In particular, the turnover of the gluon momentum distribution at = 5.0 GeV^ occurs 
below our cut-off at a; = 10~^. Hence the increase observed at = 2.5 GeV^ is no longer 
present. Instead, we see a decrease induced by the increasing value of the gluon momentum 
distribution. The numerator is approximately constant in the region of small X2 due to the 
smearing of the convolution. Since our ratio actually depends on the combination X2 + xl/2 
rather than X2, the small-X2 behavior is stabilized by taking xl finite. As displayed in Figure 
1^, our ratio saturates as X2 is reduced when is finite. In addition, its behavior is far more 
moderate. 

Our numerical analysis has shown that in a very general region the 'naive' expectation. 
Eg. (^5]) , for our matrix element is quite a good approximation. As long as xa_ and xl/xa are 
small enough, X2 + xl/2 is not too small, and 0.2 ^ Xi ^ 0.7, corrections run in the 10-15% 
level or less. We can get even better approximations by including an x^i-dependent constant 
of order 0.8. However, we must be very careful when using this approximation generally. 
It is easy to see from the above plots that the ratio varies quite quickly as one leaves the 
region of validity. 
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Q Dependence for Xj^=x^ 



Q-=2.5 GeV ' 
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FIG. 7. When xl ^ 0, the wild dependence of the ratio on X2 is suppressed as X2 decreases. 
This leads to much more moderate behavior over the whole range of X2- 

V. PHENOMENOLOGY 



In previous studies one has assumed a Gaussian spatial nuclear distribution 

exp[— (y~/2i?^)^], which leads to a phenomeno logical form for the two-parton correlation in 
nuclei, 



(52) 



Compared with the recent HERMES experimental data on nuclear modification of quark 
fragmentation function in DIS, one has extracted the constant Ca^ ~ 0.00065 GeV^. 
The strong coupling constant as should be evaluated at a scale ~ 2.8 GeV^. Such a value 
is also consistent with what is extracted from the nuclear transverse momentum broadening 
C7a2 = 0.00086 GeV^. 

Assuming now that Kq{xi,X2,xl) is a good approximation of K{xi,X2,xl), we can 
express the two-parton correlation in nuclei of Eq. (^) as 



Tqg{x,XT,XL) 



32Rt 



■{RlAl)[Q{x + xl)G{xt) + Q{x)G{xT + xl) 



-2Q{x + xl/2)G{xt + XL)e-^' 



(53) 



Here we assume a Gaussian distribution for the two-nucleon correlation function r2- For 
moderate values of a; ^ a;/, ^ a;^ « 1, one can assume Q{x + xl) ~ Q{x), the above 
expression can be approximated as 
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T^'1{x,xt,xl) = ^{R\Al)Q{x) [G{xt) + G{xt + xl)] 



2G{xt + xl/2) „_4/44 



G{xt) + G{xt + XL 



-e 



(54) 



The above derived factorization form is very close to the phenomenological one in Eq.(^), 
especially for xl <^ xt- For large xl ^ xt, when G{xt) ^ G{xt + xl), the coefficient in 
front of the exponential factor will have additional suppression as compared to the phe- 
nomenological model. However, for not so large and xl ^ xa, G{xt) ~ G{xt + xl)- In 
this kinematic region, one can then relate the parameters in the phenomenological form to 
our result within the convolution model 

CMRa^^{R\^1)G{xt) . (55) 

Here, we have reduced the nuclear radius by 1/2 in Eq. (|54D in order to match the phenomeno- 
logical form in Eq. (|52|) . With the value of Gal ~ 0.00065 from HERMES experiment and 
G{xt) ^ 3 at XT = XB{kl)/Q'^ ~ 0.01 {xb = 0.124, = 2.8 GeV^, {k^) ~ 0.25 GeV^), we 
have 

{R\Al) ^ 1.65, (56) 

which is on the order of 1 and independent of the A as expected. One can consider this 
as a qualitative agreement between our calculated two-parton correlation in nuclei and the 
experimental measurements. 



VI. CONCLUSION AND DISCUSSION 

In this paper, we have studied the generalized twist-four nuclear parton matrix elements 
K{xi,X2, Xl), both direct and momentum-crossed ones, in the framework of a Fock hadronic 
state expansion of the nuclear states. These twist-four nuclear parton matrix elements 
determine the effect of multiple parton scattering in hard processes involving nuclei, e.g., 
the nuclear modification of the fragmentation functions. Assuming that the contributions 
of higher Fock states induced by nucleonic interaction are small, we have shown that the 
leading contribution to the twist-four nuclear parton matrix elements can be expressed as 
a convolution of twist-two nucleonic off-forward parton distributions and the two-nucleon 
correlation function inside a nucleus. In the limit of extremely large nuclei, A — > oo, or 
Xa = 1/MRa — >■ 0, we have also shown that the twist-four nuclear parton matrix elements 
can factorized into the product of twist-two nucleonic parton distributions. However, we 
demonstrated that the matrix elements K are not analytic in xa at x^i = (Corrections 
around xa = cannot be expanded as powers of xa)- 

To verify the factorization approximation, we have evaluated the twist-four nuclear ma- 
trix elements numerically as the convolution of twist-two nucleonic parton distributions and 
two-nucleon correlation functions inside a nucleus, assumed to have a simple Gaussian form. 
For Xl ^ Xa and moderate xi, we found that the factorization is a good approximation 
within ~20% for large nuclei. However, the deviations become very significant for small 
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xi or X2 + xl/2, large xa, or xl ~ xa- The corrections at small X2 + are particu- 
larly large when the gluon momentum distribution G{x2 + Xl/'^) is large. Furthermore, for 
3^1,2 ^ one can no longer express the nuclear matrix elements as the convolution of 

twist-two nucleonic parton distributions. In this region, off-forward nucleonic parton dis- 
tributions, which are not known experimentally, are needed for the numerical evaluation. 
Therefore, one could conceivably use the measured nuclear effects in this kinematic region 
to constrain the nucleonic OFPD's. 

Another important nuclear effect on the parton matrix elements that we have not con- 
sidered so far is similar to the nuclear shadowing of the parton distributions or depletion of 
the effective parton distribution per nucleon. This effect is known experimentally to be large 



for small x^, at least for the quark distributions [^]. One can understand nuclear shadow- 
ing as the consequence of coherent initial scattering processes involving multiple nucleons 



T5| , p!6[| when one calculates the nuclear parton distributions in terms of nucleonic parton dis- 



tributions. In the framework of our lowest Fock state expansion, these multiple scattering 
processes will involve direct multi-nucleonic correlations. Therefore, such multiple scatter- 
ing effects in principle are higher-twist contributions and should be suppressed at very large 
Q^. In our calculation of the twist-four nuclear parton matrix elements, such higher-twist 
contributions from multiple scattering processes involving more than two nucleons should 
also contribute, but they are suppressed at large . 
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